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We review the rich phenomena associated with neutrino flavor transformation in the presence 
of neutrino self-coupling. Our exposition centers on three collective neutrino oscillation scenarios: 
a simple bipolar neutrino system that initially consists of mono-energetic and v^; a homoge- 
neous and isotropic neutrino gas with multiple neutrino/antineutrino species and continuous energy 
spectra; and a generic neutrino gas in an anisotropic environment. We use each of these scenarios 
to illustrate key facets of collective neutrino oscillations. We discuss the implications of collective 
neutrino flavor oscillations for core collapse supernova physics and for the prospects of obtaining 
fundamental neutrino properties, e.g., the neutrino mass hierarchy and O13 from a future observed 
supernova neutrino signal. 

I. INTRODUCTION 
A. Neutrino mixing and astrophysics 

Neutrinos and the phenomena associated with neutrino flavor transformation stand at the nexus of two exciting 
recent developments: the success of experimental neutrino physics; and the tremendous growth of astronomy and 
astrophysics. The former enterprise has given key insights into neutrino mass and vacuum mixing and promises more 
[1], while the latter is providing fundamental cosmological parameters and is revealing how structure and elemental 
abundances emerge and evolve in the universe [2-5]. Moreover, there is feedback between these subjects. For example, 
observations of large scale structure and the cosmic microwave background radiation currently provide our best limits 
on the neutrino rest masses [6, 7]. In fact, both the early universe and the massive star core collapse and supernova 
explosion environments can be dominated by neutrinos. Neutrino flavor transformation in each of these environments 
may give insights into astrophysics and even possibly into fundamental neutrino properties. Obtaining these insights 
will require confident modeling of neutrino flavor evolution in environments where neutrino-neutrino interactions 
produce vexing nonlinearity. In any case, the neutrino mass and mixing data already gathered from the experiments 
make a compelling case that we must solve this problem. 

Experiments and observations to date have established that the neutrino energy (mass) states {vi) {i = 1,2,3) 
are not coincident with the weak interaction states {va) (ct = e, M,t). The relation between these bases is given 
by |i/q.) — J2i ^ai Wi) [^]' where the Maki-Nakagawa-Sakata (MNS) matrix elements Uai are parameterized by three 
vacuum mixing angles {9i2, O23, ^13) and a CP- violating phase S. Two of these are measured outright, sin^ 612 « 0.31 
and sin^ 6*23 ~ 0.5, while there is a firm upper limit on a third, sin^ 613 < 0.04 at 2a [9]. Observations of solar neutrinos 
show flavor conversion in the ^ ^m/^ channel with a characteristic mass-squared splitting Atoq « 7.6 x 10~^eV^. 
Atmospheric neutrino measurements show near maximal vacuum mixing in the i^^ channel with corresponding 

mass-squared splitting Am^tm ~ 2.4 x 10"'^ eV^. However, experiments do not reveal the absolute neutrino rest 
masses nii (i = 1,2,3), nor do they show whether these neutrino mass eigenvalues are ordered in the normal mass 
hierarchy (7713 > TO2 > mi) or in the inverted mass hierarchy (7712 > nii > ms). 

Future terrestrial neutrino experiments [10-15] will be directed primarily towards measuring ^13 and, if ^13 is big 
enough, the neutrino mass hierarchy and possibly even the CP- violating phase 6. The planned reactor and long 
baseline experiments may be able to measure ^13 if it satisfies sin^ ^13 > 10"'^. This limit is set ultimately by 
constraints on the neutrino flux and detector mass. To find significantly larger neutrino fluxes we must turn to cosmic 
sources, e.g., core collapse supernovae [16-23]. 

Stars more massive than ~ 8 Mq end their lives in gravitational collapse and the production of a neutron star 
or, if they are massive enough, a black hole remnant [24-27]. Neutrinos play a role in nearly every aspect of the 
evolution of these core collapse supernovae, from dominating lepton number and entropy loss from the epoch of core 
carbon/oxygen burning onward, to providing the bulk of energy and lepton number transport during collapse itself, 
and even to providing the heating necessary to engender convection [28-33] and, e.g., the Standing Accretion Shock 
Instability (SASI) which may create an explosion [34, 35]. 

A key point is that gravitational collapse causes an appreciable fraction of the rest mass of the Chandrasekhar-mass 
1.4 M0) core to appear as seas of trapped neutrinos which subsequently diffuse out of the core on time scales of 
seconds [36-38]. At core bounce, the energy in the neutrino seas trapped in the core is ^ 10*^2 ergs, but by - 10 s 
after core bounce some 10^'^ ergs, or ~ 10% of the rest mass of the core has been radiated away as neutrinos of all 
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kinds. The emergent neutrino and antineutrino energy spectra and fluxes at the neutrino sphere vary with time 
post-core-bounce, but there will be epochs where the energy spectra and/or the luminosities in the various neutrino 
flavors will differ. Moreover, the charged-current, flavor specific neutrino interaction processes such as j/g ^ p + e^ 
and + P ^ n + are important both for energy and electron lepton number deposition as well as for determining 
neutrino transport physics [39-41]. It is therefore interesting and necessary to assess whether interconversion of 
neutrino flavors in the supernova environment affects explosion physics and neutrino-heated nucleosynthesis and how 
such flavor transformation might affect a supernova neutrino burst signature in a terrestrial detector. 

B. Brief history of collective neutrino oscillations 

Early studies of neutrino flavor transformation centered on solar neutrinos, especially after it was recognized flrst by 
Wolfenstein [42] , and then Mikheyev and Smirnov [43] , that the medium through which the neutrino propagates could 
alter the effective neutrino mass and mixing properties. Almost immediately after the Mikhyev-Smirnov- Wolfenstein 
(MSW) mechanism was understood, it was realized that the forward coherent scattering of neutrinos with other 
neutrinos, somewhat misleadingly called "neutrino self-coupling" or "neutrino self-interaction" , could generate a 
similar effect [44, 45]. Subsequently, the effects of neutrino self- interaction were investigated independently in the 
core collapse supernova and the early universe scenarios. The studies of the supernova environment flrst focused on 
MSW- like effects [39, 46-49]. However, it was pointed out that the neutrino self- interaction potential is very different 
from the matter potential in that it can have non- vanishing off-diagonal elements in the interaction basis [50, 51]. 
Studies using the complete neutrino self-interaction potential showed that neutrinos could experience "self-maintained 
coherent oscillations" or "collective oscillations" in lepton-degenerate early universe scenarios [52-62]. By collective 
oscillations we mean a significant fraction of neutrinos oscillate coherently with respect to each other. 

Early research on MSW-like evolution, where the flavor off-diagonal potentials were minimal, as well as the flrst 
paper [63] to point out that collective effects could occur in supernovae all focused on relatively large neutrino 
mass-squared differences. However, it was eventually realized that neutrinos and antineutrinos could be transformed 
collectively and simultaneously even for the small, measured neutrino mass-squared differences [64, 65]. Additionally, 
Reference [66] showed that ordinary matter does not necessarily suppress collective neutrino oscillations, at least in 
homogeneous and isotropic environments. It was then demonstrated that collective neutrino oscillations indeed can 
occur in a spherically symmetric supernova model [67, 68]. These works showed that the neutrino energy spectra 
would be modifled differently for the normal and inverted neutrino mass hierarchies (see Figure 1). Over the 
past few years, many papers have been written on the collective neutrino oscillation phenomenon and its physical 
effects [69-108]. Following this literature can be bewildering. For example, the synchronized and bipolar neutrino 
oscillations are frequently perceived as the two most important collective neutrino oscillation modes in supernovae. 
However, numerical studies suggest that supernova neutrinos are probably never synchronized because of the non- 
vanishing matter density [74], and the bipolar oscillation is not a collective neutrino oscillation mode in anisotropic 
environments such as supernovae [70]. 

C. Goal and organization of this review 

The purpose of this article is to provide a relatively short but in-depth review of the properties of collective neutrino 
oscillations that are reasonably well understood. In particular, we wish to elucidate the physics behind the striking 
features such as the swaps of supernova neutrino energy spectra caused by collective oscillations as shown in Figure 1. 
Our strategy is to illustrate key facets of collective neutrino oscillations by describing three, increasingly more complex 
models. 

The rest of the review is organized as follows. In Section II we try to make the connection between the "wavefunction 
language" and the "spin language", commonly used for studying non-collective and collective neutrino oscillations, 
respectively. In Section HI we discuss the flavor evolution of a homogeneous and isotropic neutrino gas that consists 
initially of mono-energetic pure Ve and Ue- This simple model can be solved analytically and offers important insights 
into the behavior of more complicated systems. In Section IV we discuss homogeneous and isotropic gases that consist 
of neutrinos with continuous energy spectra. We explore adiabatic solutions to the neutrino flavor evolution equations 
and demonstrate how the spectral swap/split phenomenon can be explained by one of these solutions. In Section V 
we discuss some of the important neutrino oscillation properties that are unique to anisotropic environments and 
explain why the spectral swap/split phenomenon may occur in these environments despite the anisotropy in neutrino 
flelds. In Section VI we apply the current understanding of collective neutrino oscillation phenomenon to supernova 
environments. In Section VII we give a summary and point out several issues in collective neutrino oscillations that 
remain to be understood. 
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FIG. 1: Survival probabilities P^^ for neutrinos (left panels) and antineutrinos (right panels) as functions of both neutrino 
energy E and emission angle i9_r in a numerical calculation using the neutrino bulb model (see Figure 3) and the two-flavor 
mixing scheme. The most prominent features on this figure are the approximately angle-independent step-like changes in the 
neutrino survival probabilities (left panels) — these are the spectral swaps/splits discovered in Reference [68]. The energy 
spectra of Ve and t'^ with energy below (above) £ ~ 9 MeV in this calculation are almost completely swapped in the upper 
(lower) panels which employ the normal (inverted) neutrino mass hierarchy. Spectral swaps/splits are the result of a collective 
neutrino oscillation mode, the precession mode, which is one of the main topics of this review. The vertical fringes (the energy- 
dependent features in the figure) are the result of MSW flavor transformation which is energy dependent. The horizontal fringes 
(the angle-dependent features) are the result of the kinematic decoherence of bipolar neutrino oscillations [70] (see Section V A). 
The movie version of these calculations is available in Reference [85]. Figure adapted from Figure 3 in Reference [67]. Copyright 
2006 by the American Physical Society. 

II. NEUTRINO MIXING IN DENSE NEUTRINO GASES 
A. Equations of motion 

Here we focus on two-flavor neutrino mixing scenarios, i.e., a = e, fi and z = 1, 2, where is a linear combination 
of the physical and |i/t-). We do this for pedagogical purposes, although this is physically justifiable because the 
physical and arc nearly maximally mixed in vacuum and experience nearly identical interactions in the supernova 
environment [109, 110]. Wc will discuss collective neutrino oscillations with the full three-flavor mixing machinery in 
Section IV C. We consider only coherent neutrino flavor evolution, where the effects of neutrino inelastic scattering 



4 



and associated quantum decoherence can be neglected. This will be generally applicable in the region well above the 
neutron star in supernova models. (Solution of the complete problem of neutrino flavor evolution with both elastic 
and inelastic neutrino scattering would necessitate the use of the full quantum kinetic equations [7, 41, 111-113].) We 
also assume that neutrinos are relativistic and that general relativistic effects can be ignored. With these assumptions, 
a mean-field Schrodinger-like equation 



i— -0 = Hip 
dx 



(1) 



is taken to describe flavor evolution along a neutrino world line [114]. Here, we take h = c = 1, x is the distance along 
the world line of the neutrino, ip is the neutrino flavor wavefunction, taken to be a vector in (neutrino) flavor space, 
and 
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+ V2GFn,l + V2Gf J d3p'(l - p • p')(pp' - Pp- 



(2) 



is the Hamiltonian. (See References [115-119] for discussions of the applicability of the one-particle effective ap- 
proximation assumed in Equation (1).) The wavefunction for an antineutrino also obeys Equation (1), but with the 
replacements A — >■ —A and H^^, — > — H*j^ in Equation (2). The flavor space is spanned by the neutrino interaction 
basis or, equivalently, spanned by the mass basis. We will use the word "flavor" in a more general sense and we will 
avoid the phrase "flavor state", which in the literature can be taken to mean either the interaction state l^a) or the 
flavor quantum state lip) . (We will use symbols in the sans serif font, e.g., H, to denote basis-dependent matrices in 
flavor space.) 

The first term in Equation (2) induces neutrino flavor transformation because in the interaction basis 



B = U ( idiag[-l, 1] ] Ut 



- cos 26'v sin 20^ 
sin 26'v cos 29y 



(3) 



is non-diagonal. In Equation (3), U is the MNS matrix, and 9^ is the so-called vacuum mixing angle and is within 
the range (0, 7r/4]. In Equation (2) Am^ is the mass-squared difference appropriate for [1/2) and and E is the 
energy of the neutrino. Here Am^ > and Am^ < correspond to the normal neutrino mass hierarchy (NH) and the 
inverted neutrino mass hierarchy (IH), respectively, as discussed above for the full three-flavor mixing case. Although 
this formalism is completely general, one of the most interesting supernova cases is where |Am^| ~ Am^^j^ and 
e»v ~ 6*13 < 1. 

The second term in Equation (2) arises from coherent neutrino-electron forward exchange scattering [42]. In this 
term, which is referred to as the matter term, Gp is the Fermi constant, rie is the net electron number density, and 
L = diag[l, 0] in the interaction basis. Note that adding/subtracting a multiple of the identity matrix to/from H (e.g., 
the contribution of neutral-current neutrino-electron scattering) gives only an overall phase to ip and, therefore, does 
not affect neutrino oscillations. 

The last term in Equation (2) stems from coherent neutral-current neutrino-neutrino forward exchange scattering 
[44-46, 50], where p and p' are the unit vectors for the propagation directions of the test neutrino and the background 
neutrino or antineutrino, respectively. In the interaction basis, at location x and at time t, the (flavor) density matrices 
for neutrinos and antineutrinos with momentum p' and with our assumptions can be written as 



[pp'{t,X.)]ap = ^nj./,p.(t,x)(l/a|l/V',p'(t,x))('i/'^'^p/(i,x)|z^^), 

[pp'{t,x.)]pa = ^np/,p.(i,x)(Pa|V'p',p'(i,x))('(/;p',p'(i,x)|j?/3), 



(4a) 
(4b) 



respectively, where [V'!^'(p').p') is the state of a neutrino (antineutrino 9') with momentum p', and n^'_p' (rip'.p') 
is the corresponding number density of the neutrino (antineutrino). Note that the order of the indices on the matrix 
representation of p in Equation (4b) follows the convention in Reference [111]. The advantage of this definition is 
that p will transform in the same way as does p when transforming from the interaction basis to the mass basis or 
vice versa. 



B. Neutrino flavor polarization vector 



It becomes more difficult to analyze neutrino oscillations using the wavefunction formalism when H^^, is significant. 
This is because H^^ is a sum of density matrices which involve bilinear forms of the wavefunctions. The density 
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FIG. 2: Geometric pictures for vacuum oscillations and for MSW flavor transformation. The flavor Hilbert space spanned by 
\va) {a — e,/i) or (i — 1,2) can be mapped onto the flavor Euclidean space that is spanned by Sa' or ei^' (a = 1,2,3). 
The interaction basis vectors e^^g can be obtained by rotating the (vacuum) mass basis vectors e^^^ by 26^ about ej^' = Cj^' 
. The left panel shows the polarization vector P„, which describes a neutrino initially in and its precession about B. 

The projection of this precession motion onto the 63^' axis represents the flavor oscillation of the neutrino. The right panel 
shows the precession of in the presence of ordinary matter. If the matter density varies slowly, the angle between and 
H = ujB + XL remains constant. This represents adiabatic MSW flavor transformation of the neutrino. See also Reference [120] 
for a more detailed discussion on this geometric interpretation. 



matrices contain all the physical information for the neutrino mixing problem. The diagonal elements of the density 
matrices give the number densities of neutrinos in the weak interaction states or the mass states, depending on 
the basis used, and the off-diagonal elements of the density matrices contain the neutrino mixing information. For 
simplicity let us first take a working example, a homogeneous and isotropic neutrino gas whose flavor content can vary 
with time t. Homogeneity and isotropy implies that the factor (1 — p-p') in Equation (2) averages to 1. Therefore, the 
neutrino propagation direction does not matter, and pp{t,x.) — > pE{t) and pp(t,x) — > pE{t). These density matrices 
obey the equations of motion (EoM): 
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B + AL + V2Gf / dE'ipE' - pe'), Pe 
Jo 

POO 

B + XI + V2Gf / dE'ipE' - PE'), PE 
Jo 



(5a) 
(5b) 



Because pE and pE are 2x2 Hermitian matrices, they can be mapped into vectors in a three-dimensional Euclidean 
space which we will also call flavor space. We define the components of the (neutrino flavor) polarization vector P„ 
to be 

J_\ / |Ato^| \ ^ r Tr(pB CTa) for neutrino, 
n,yj\ 2a;2 / \~Tr{pEcra) for antincutrino, 

where the (angular) vacuum oscillation frequency is a; = -^g- for neutrinos and (^f^ for antineutrinos, and CTq 

(a = 1,2,3) are the Pauli matrices. (We use symbols with the vector hat, e.g., P, to denote a vector in flavor space 
and symbols in the bold font, e.g., p, to denote a vector in physical three-dimensional coordinate space.) We note 
that Puj can be normalized by an arbitrary factor. For example, if Pj is normalized to unity, then pE oc {1 + ^Puj ■ 
for neutrinos and ex (1 — ^Pui ■ f?) for antineutrinos. In Equation (6) we defined Pj to be normalized by n^, 
the initial total number density of a certain neutrino species 1^ which is chosen to be Dg in the rest of the review. 
However, we usually take p = V^Gpn^ as a measure of the strength of neutrino self-interaction. Therefore, in some 
cases it is more appropriate to normalize P^i by the number density of other neutrino species, e.g., when rip^ is 
negligible. If neutrinos and antineutrinos are all in the interaction states, then P^ oc (^^jy^.o; — 't-j^^ ,1^)63^^ for neutrinos 
and Pj cx — (rip^.tj — np^_i^)e'^^ for antineutrinos, where n^.^ {ly — iye,i^p,,i^e,i^p,) is the corresponding number density 

in the neutrino or antincutrino mode w, and 63^'' is one of the interaction basis vectors. 

Since the Pauli matrices are traceless, the trace of the density matrix is not contained in the polarization vector. 
According to Equation (5) the traces of the density matrices do not change with time. (This corresponds to one of 
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our assumptions that neutrinos are not created or annihilated.) These terms can be easily reintroduced, e.g., for the 
calculation of the neutrino energy spectra. Using Equations (5) and (6) it can be shown that 

= {ujB + \L + ^iD) X P^, (7) 

where vectors B = Tr(B(7) and L = Tr(LtT) are parallel to the (vacuum) mass and interaction basis vectors 63^'' 
and respectively (see Figure 2), and D = Pujduj is the total polarization vector. In the absence of neutrino 
self-coupling, P^ can be thought of as a "magnetic spin". In this analogy the "magnetic spin" is coupled to two 
"magnetic fields", B and L, with gyromagnetic ratios — w and —A, respectively. (Note that a real magnetic spin s 
with the gyromagnetic ratio 7 in the presence of magnetic field B obeys EoM s — — 7B x s.) Equivalently, behaves 
like a "magnetic spin" coupled to the total "magnetic field" H = ujB + XL with "gyromagnetic ratio" —I. This picture 
allows geometric interpretations for both vacuum oscillations and MSW flavor transformation (see Figure 2). 

A few comments about the polarization vector notation are in order before we consider neutrino self-coupling, (a) 
The polarization vector notation is fully equivalent to the neutrino flavor isospin (NFIS) notation [66] where the NFIS 

for a neutrino or antineutrino is s^j = ^7^7. (b) In Equation (6) is defined with a minus sign for the antineutrino. 

Although the physics is not changed by the choice of notation, our definition of P^ is convenient in analyzing collective 
neutrino oscillations when both neutrinos and antineutrinos are present (see Section II C). This notation has been 
adopted in most recent literature. 



C. Synchronized neutrino oscillations, corotating frames and matter effects 

The neutrino self-coupling is represented as the coupling between polarization vectors in Equation (7). To have 
a feeling for the effect of this coupling, let us consider a homogeneous and isotropic neutrino gas with A = and 
/i = const. Using Equation (7) it can be shown that the "total energy of the magnetic spins" 

£^ f ujiP^- B)dLu + ^D^ (8) 

is constant in this case [66], where the first term is the "total energy" of coupling between the "magnetic field" 
and "spins", and the second term is the "total spin-spin couphng energy". If /i is large, then ~ yj2£ j [i is 
approximately constant. This implies that a dense neutrino gas can experience "self-maintained coherent oscillations" 
[56]. For example, if a dense neutrino gas consists initially of neutrinos of the same flavor (so that all P,^ are initially 
aligned), then the flavor evolution of these neutrinos is coherent (i.e., all P^ remain aligned) even when they have 
different energies. This phenomenon has been termed synchronized neutrino oscillations, because all neutrinos (and 
antineutrinos) in such a system oscillate collectively with (angular) frequency fisync- The synchronized oscillation 
frequency is an average of all w's [60]: 

/oo 
[p ■ P^)ujduj. (9) 
-00 

Consider another neutrino system that is similar to the synchronized system discussed above except that the 
oscillation frequency lo of each neutrino or antineutrino is shifted by a common value wq. Any polarization vector, say 
P^j, in this system should move in a way similar to Pu)-^)^ the synchronized system, except for an extra precession 
about B with frequency wq. In other words, this system behaves just like the synchronized system in a reference 
frame that rotates about B with frequency wq- Indeed, in this non-inertial corotating frame [66], each polarization 
vector is coupled to a non-physical held —ujoB and Equation (7) (with A = 0) becomes 

P^^[iio-ujo)B + fiD]xP^. (10) 

Therefore, this neutrino system experiences synchronized flavor transformation, just like the one discussed above, but 
with rJgync shifted by luq. Because the oscillation frequency lu can be shifted to any value by using an appropriate 
corotating frame, the criterion for synchronization should not be ^ 3> but rather ^ ^ Alo [66], where (uj) and 
Auj are the average value and the spread in uj for the neutrino system, respectively. Note that the reversal of the 
direction of Pi^ for the antineutrino (Equation (6)) is important when the corotating frame is used. If this reversal 
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FIG. 3: The geometric layout of the neutrino bulb model. In this model all neutrinos are emitted half-isotropically from the 
surface (neutrino sphere) of the PNS which has radius R. Spherical symmetry and isotropic emission on the neutrino sphere 
imply that all neutrinos with the same initial flavor, energy and emission angle have identical flavor evolution histories. 
In this model the neutrino polarization vector P[^,,j(r) is uniquely determined by uj, 1} (or and r. Here -d is the angle at 
radius r between the neutrino trajectory direction and the radial direction. Figure adapted from Figure 1 in Reference [68]. 
Copyright 2006 by the American Physical Society. 



is not incorporated into the definition of P^j, then the direction of also changes if the sign of the corresponding 
vacuum oscillation frequency changes on transformation to the corotating frame. 

Note that if B is parallel to L, the matter effect can be completely removed by transforming to an appropriate 
corotating frame. For a general case, on transforming to the corotating frame Equation (7) becomes 

P^ = {ujB + i^iD) X P„ and B = -XL x B. (11) 

The matter effect does not disappear here, but rather causes B to rotate in the corotating frame. If the matter density 
is large (A » |w|), however, the fast-rotating B in Equation (11) can be replaced by {B ■ L)L for collective neutrino 
oscillations [66]. In other words, for collective neutrino oscillations and in the presence of large matter density, matter 
effects may be "ignored" and the effective neutrino mixing parameters become 9f,a — and Am^g = Am^ cos20v 



D. Solving for supernova neutrino flavor evolution 

Although collective neutrino oscillations may occur in any environment where neutrino fluxes are significant, recent 
studies of this phenomenon have focused on the core collapse supernova environment. The supernova environment 
is far more complex than the early universe, in part because of its inhomogeneity and anisotropy. This complexity 
is enhanced because neutrinos of different flavors and energies and propagating in different directions arc coupled 
by neutrino self- interaction [48]. Full simulations of neutrino oscillations with neutrino self-interaction in a general 
supernova environment are beyond current numerical capabilities. Here we briefly discuss two schemes commonly 
used in investigating collective neutrino oscillations. Both of these schemes employ the neutrino bulb model where 
the supernova environment is spherically symmetric around the center of the proto-ncutron star (PNS) (see Figure 3). 
The polarization vectors in this model obey EoM 

cosi?-^/l,^(r) = + A(r)L + i?^^,^(r)] x (12) 
dr 

where 

/oo />1 
dw' / d(cosz?')(l-cosz?cosi?')^l',<?'W- (13) 

Note that the maximum value of 'dir) is i^max — arcsin(i?/r) in the neutrino bulb model (see Figure 3). Also note 
that here we choose to normalize P^ by rivS^) = Lp^/{2TrR'^{Ep^)), the total number density of Df. at the neutrino 
sphere, where Lp^ and {E^J are the energy luminosity and the spectrum-averaged energy of Pe at the neutrino sphere, 
respectively. For the neutronization-burst epoch where the i^e flux is much larger than the fluxes of all other neutrino 
species, P^j should be normalized by n^^ (R) instead. In typical numerical simulations, R is in the range of 10-60 
km, and the luminosities and the average energies of neutrinos are in the ranges of 10^°-10^"^ erg/s and 10-30 MeV, 
respectively. 

Equation (12) can be solved numerically without any further assumptions. This is the "multi-angle scheme". The 
other scheme is the so-called "single-angle scheme". In this latter scheme it is assumed that Pu:,^{r) = Pi^{r) is the 
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same for different neutrino trajectories. There are several variants of the single-angle scheme. These variants lead 
to qualitatively similar results. In one of the variants Pu:{r) is computed along the radial direction {'dn = 0) and 
Equation (12) becomes 



d 
dr 



P^{r) ^ [loB + \{r)L + ^/2Gvni,V{r / R)D] x it(r), 



(14) 



where the geometric factor 'D{r/R) 



1 — A /I — (■^)^ partially accounts for the angle effect and geometric 



dilution of the neutrino fluxes in the neutrino bulb model. Comparing Equations (7) and (14), it is clear that the 
flavor evolution of neutrinos in the single-angle scheme is equivalent to that of a homogeneous and isotropic neutrino 
gas expanding with "time" r. In this analogy the strength of the neutrino self-coupling is ^(r) = ^/2GYnv^{R)'D{r / R). 
The radial direction is a rather special direction in the neutrino bulb model. In another variant of the single-angle 
scheme it is assumed that all neutrinos are emitted with 'd^ = 7r/4 [73]. Alternatively, H,jy^^ can be averaged over 
neutrino trajectories [68, 88]. Each of these variants also leads to Equation (14) when r ^ i?. However, in this limit 

'D{r/R) — s- i (-^) for these variants instead of ^{^)^ in the single-angle scheme which employs the radial trajectory. 
These variants can improve the agreement between the results of the single-angle and multi-angle calculations [88] . 

Although the single-angle scheme is frequently used for its simplicity, it must be emphasized that it misses some 
of the important properties (e.g., anisotropy) of the neutrino bulb model and can lead to incorrect results when such 
properties play important roles in collective neutrino oscillations (see Section V). 



III. SIMPLE BIPOLAR NEUTRINO SYSTEMS 

A. Bipolar systems and the flavor pendulum 

To illustrate another important example of collective neutrino oscillations, let us consider a homogeneous and 
isotropic gas that initially consists of mono-cnergetic Ve and Pe ■ This neutrino system is represented by two polarization 
vectors, Pi^ and P-i^ for the neutrino and the antincutrino, respectively. We assume that A = 0, /i is fixed, and that 
l-Pc^l = (l+e)|-P-Lj|, where e is the fractional excess of neutrinos over antineutrinos. At i = 0, P^ points in the direction 
of e'^\ which is tilted away from 63^^ by 29^, and P^i^ points in the direction opposite to P^. Neutrino systems that 
are represented by two nearly oppositely directed polarization vector groups are called "bipolar systems" . Using the 
corotating frame technique, the discussion in this section easily can be applied to, e.g., a gas consisting initially of Ve 
and Vfj^ with energies E^^ ^ Ei,^ . 

A peculiar case is where e = and initially is aligned with eg^-* = —B (i.e., 9^ = 0). Using energy conservation 
(Equation (8)) it can be shown [66] that if a; > 0, neither P^ nor P^^i will move, and so the initial configuration of the 
system is absolutely stable. On the other hand, if w < and /j, ^ \uj\, P^j and P-i^ can nearly swap their directions 
(but with a slight bend towards each other). This implies that the initial configuration of the system is unstable. 
Therefore, the v^-Ve, system may experience insignificant flavor oscillations when Am^ > and ^ 1. However, 
this system can experience significant flavor oscillations when Am^ < and 9^ <^\. Collective neutrino oscillations 
of this kind are known as "bipolar oscillations" . 

The EoM of the simple bipolar system 

it = (wP + [iD) X /L, and = {-ujB + ^D) x (15) 

have been solved analytically [56, 57]. Instead of presenting this solution, let us rewrite [69] Equation (15) as 

D = M^^g X g and D = i-L^^q x q + asQ, (16) 

where q = Q/\Q\ = (P^ — P-uj ~ ^^)/IQl! 9 = ^f^^\Q\B^ and = q ■ D is constant. Equation (16) describes the 

motion of a flctitious gyroscopic pendulum, or "flavor pendulum" , with total angular momentum D in a uniform 
gravitational fleld where the acceleration of gravity is g. The pendulum consists of a massless rod with a point 
particle of mass /z"^ and spin CTs attached to the end of the rod at position q. We note that the flavor pendulum that 
represents a symmetric bipolar system (e = 0) has no internal spin. The stable and unstable configurations of the 
system discussed above correspond to the lowest and highest positions, respectively, that the pendulum can reach. 

Generally, the flavor pendulum can experience two kinds of motion: a precession about the B axis and a nutation 
around the average precession track. The nutation motion corresponds to bipolar neutrino oscillations. However, like 
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FIG. 4: Evolution with decreasing for a simple bipolar system that consists initially of mono-energetic and v^. In this 
example the neutrino excess is £ = 0.25 and the mixing parameters are — 0.01 and Am^ < 0. The left panel shows the 
configuration of the polarization vectors, where Qv {9i?) is the angle between _fL (P-i^) and B = —63^'. The thick lines are 
for an assumed pure precession motion of the flavor pendulum, and the thin lines are for a case where fi decreases slowly and 
linearly with time. At /x > /icr the flavor pendulum is a "sleeping top" and there are no flavor oscillations. When fi < 
the flavor pendulum experiences precession as well as nutation around the average precession track. The right panel shows 
the precession frequency Q, of the flavor pendulum assuming that it is executing pure precession. In this case, the value Q 
approaches f2sync and u in the limits /j. ^ 00 and /i — >■ 0, respectively. The analytical solution for the pure precession motion 
of the polarization vectors can be found in Reference [77]. This flgure is adapted from Figures 3(a) and 4 in Reference [71]. 
Copyright 2007 by the American Physical Society. 



a child's top, the flavor pendulum can "defy gravity" and precess almost uniformly if ^ is large enough. In particular, 
the flavor pendulum can become a "sleeping top" and will not fall from its highest position if [69, 71] 

(vi + e - ly 

This precession behavior of the flavor pendulum in the large limit represents synchronized oscillations of the bipolar 
system with ilsync = (1 + 2e^^)Lo. 



B. Bipolar systems with slowly decreasing neutrino density 

Let us now focus on the IH case with 6^ <^ 1. Signiflcant neutrino oscillations can occur in this case. If the neutrino 
density decreases, the mass fi~^ of the pendulum becomes larger. The swing amplitude of the pendulum that represents 
a symmetric bipolar system decreases as fi decreases. The maximum swing amplitude of this pendulum can be found in 
the adiabatic limit where ^ changes slowly [71]. In this case Puj and P^i^ become aligned and antialigned, respectively, 
with B (i.e., Ve — ^ 1^2 and i>e — )■ P2) as fj, decreases toward 0. The asymmetric bipolar system is more interesting. 
Neutrino oscillations in this case are synchronized and the flavor pendulum precesses uniformly in the limit n ^ 
If jU decreases very slowly, the flavor pendulum can still experience a nearly pure precession motion for any given 
11. In this case Pui and P-ui He in the same plane with B and their directions can be readily solved for [71, 77] (see 
Figure 4). We note that, assuming e > 0, P-oj becomes antialigned with B (i.e., De V2) as — >■ 0. Meanwhile the 
direction of P^, can be determined from the constancy oi D ■ B [69] (see Equation (16), noting that g B). 



IV. ISOTROPIC AND HOMOGENEOUS NEUTRINO GASES 
A. Static solutions for the neutrino flavor evolution equation 



We next consider a homogeneous and isotropic gas that consists of neutrinos with continuous energy spectra. Let 
us flrst seek static solutions to the EoM for the polarization vectors (Equation (7)) in the case where both the matter 
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density and the neutrino densities are constant. One possibility is that the solution is stationary and the P^i's do 
not evolve with time t. This is only possible if the "spin" is parallel to the corresponding total "magnetic field" 
= ujB + XL + fiD. For this case, and noting that does not depend on P^, we can obtain the following equations 
for the total polarization vector D [74] : 



Di = (Asin26'v + fiDi) f ^^^^dw, (18a) 



^Ld I PuJ I 



D2 = fiD2 I ^^V^dw, (18b) 



D3= (-cj + Acos26>v + A*£'3) L ' duj, (18c) 

J — oo \-^ui\ 

where Da — D ■ ei^^ (a = 1,2,3), and = +1 (—1) if Puj is aligned (antialigned) with i?^. The stationary solution 
to the EoM for P^^ can be found from Equation (18) and the alignment condition. We note that Equations (18a) and 
(18b) generally imply that D2 = if A 7^ which, in turn, implies that • 62^'' = for any uj. 

Equations (18a) and (18b) become equivalent when A = 0, and D is undcrconstraincd by Equation (18) in this case. 
This is because when A = 0, the EoM for Puj possesses a rotational symmetry about B. In other words, if {P(^(i)| Vw} 
solves the EoM, then {P^(t)| Vw} also solves the EoM, where P[j{t) is obtained from Puj{t) by rotation about B by an 
arbitrary angle (j). Here (j) is independent of to and t. This is a generalization of the rotational symmetry of the flavor 
pendulum about the "gravity" vector g oc B. Such a symmetry generally implies the existence of a collective motion 
and a conservation law. For example, a translational symmetry of a group of particles along some direction implies 
the possibility of the collective motion of the particles in that direction and the conservation of the total momentum 
of the particles in the same direction. In Section III it was shown that this rotational symmetry implies the possibility 
of a pure precession of the flavor pendulum in which both polarization vectors of the simple bipolar system process 
with the same frequency fi. Therefore, it is natural to seek a static solution to the EoM with A = 0, and fi constant, 
and in which all the P^j's precess about B with the same frequency n. In this case all the Pi^'s are stationary in a 
corotating frame which rotates about B with frequency f2. In this corotating frame, P^j is either aligned or antialigned 

with = (cj — Q)B + fxD, and the components of D in the corotating frame can be found from equations similar to 
Equation (18). These equations can be recast as the following two simple sum rules [72]: 

1 = r , dc., (19) 

n^r ^ ---'^"1 d., (20) 

where D± is the component of D that is perpendicular to 63^^ The rotational symmetry of the system about B 
implies that = —D ■ B is constant. 



B. Adiabatic solutions and the spectral swap/split 



If A and fi vary slowly with time, adiabatic solutions can be obtained which correspond to the static solutions 
discussed above. We shall term these the (adiabatic, ;/-enhanced) MSW solution and the (adiabatic) precession 
solution. Figure 5 shows these two adiabatic solutions for a single-angle scheme together with the corresponding 
numerical solutions for Equation (14). In these calculations it is assumed that the most abundant neutrino species at 
the neutrino sphere are and Pe, and that 9^ ^ 1- For the precession solution the matter field is "removed" using 
an appropriate corotating frame. The fiavor pendulum model provides insight into how the numerical solution evolves 
from the MSW solution towards the precession solution. In the IH case a large matter density essentially keeps the 
flavor pendulum (with e > 0) near its highest position. This configuration becomes unstable when fi < /Xcr- In the 
NH case the z^-enhanced MSW flavor transformation has the effect of raising the flavor pendulum to near its highest 
position, and the configuration is again unstable when fi < 

Just as in the conventional adiabatic MSW fiavor transformation case (Figure 2), in the adiabatic precession 

solution P^j follows Hui which changes its direction (and magnitude) as decreases. This induces neutrino fiavor 
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FIG. 5: Comparison of the numerical solution, the //-enhanced MSW solution, and the precession solution in a single-angle 
scheme, where 6^, = 0.1, (sx) oc D± and (sz) oc D3. The numerical solution follows the MSW solution at first. In the IH case 
the MSW solution becomes unstable at r ~ 88 km in this calculation and thereafter the numerical solution shows oscillations 
around the precession solution. In the NH case the numerical solution follows the MSW solution through the resonance (where 
{sz) ~ 0) before it shifts to follow the precession solution track. Figure adapted from Figure 1 in Reference [74]. Copyright 
2007 by the American Physical Society. 



transformation. In particular, as /i — > 0, H^^ — > (w — fio)^, where f7o = ^^(m = 0)- This means that the adiabatic 
collective precession mode converts the initial Ve into the mass state jz^i) or \v2) depending on whether w is smaller 
or larger than J7o [68]. This phenomenon, the "stepwise spectral swap" or "spectral split", is most dramatic when 
6'v ^ 1 (see Figure 6). The swap/split energy Eg = \ | can be determined from the constancy oi D ■ B [72]. 



C. Precession solution in the three-flavor mixing scenario 



The neutrino polarization vector defined in Equation (6) can be generalized easily to the three-flavor mixing scenario 
by replacing the Pauli matrices with the Gell-Mann matrices (a = 1, 2, . . . , 8) [81, 120]. However, because an eight- 
dimensional polarization vector, or Bloch vector, cannot be as easily visualized as its three-dimensional counterpart, 
we will discuss the collective precession mode by using the matrix formalism. To this end, we define the polarization 
matrix P^^ = \ Y^^a=i^-^^ 'where P^j.a is the ath component of the Bloch vector P^^. The polarization matrix 
obeys the EoM 



iP,^ = [wlBl+whBh + ^D, P,^], 



(21) 
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FIG. 6: Illustration of the "stepwise spectral swap" phenomenon in the two-flavor mixing case with Am^ < (IH) that was 
discovered in Reference [68]. The left panel shows the stepwise swapping of Ue. and energy spectra about iJs — 9 MeV 
in a single-angle scheme. The spectra of Ve and are nearly fully swapped in this calculation. The right panel shows the 
corresponding survival probability P^u which is a step- like function of uj. For the NH case the step-like structure of Puvij-o) 

is pushed rightward to fio > 0. Because Es = 
phenomenon is also sometimes called "spectral split 



I splits a neutrino spectrum into two parts with different flavors, this 
Reprinted with permission from H. Duan, AIP Conf. Proc. Vol. 1182, 



Page 37, 2009 (Reference [121]). Copyright 2009, American Institute of Physics. 



where D = Pc^do; is the total polarization matrix. In Equation (21), ojl = and Bl = — 5A3 (in the mass 

basis) correspond to the small mass splitting which we define as 5m? = m\ — m\ ~ Amg. Also in Equation (21), 
wh = o; = i'^^ and Bh = — ^Ag (in the mass basis) correspond to the large mass splitting which wc define as 

A7Ti^ = m§ — \{'^\ + "^2) — ^^'^atm- -f^^r simplicity we have ignored the matter field since it can be removed by 
employing the corotating frame technique [82]. 

The static precession solution (with constant /i) can be obtained by assuming that {P^n f^Lj ^^h| Vw} solves equation 

[H^, Pu] = [K - f^L)BL + (c^H - f^H)BH + mD, p.] = 0, (22) 

where P^^, f^L and JIh arc constant, and D = P^dw. In the two-flavor mixing scenario (Bl = 0), Equation (22) 
corresponds to the condition that the "spin" P^^ is stationary in a corotating frame and is parallel to the total 

"magnetic field" Hui in this reference frame. Equation (22) is called the precession ansatz [82] because it implies a 
static precession solution to the EoM (Equation (21)). This solution, {Ptj(i), f^L, ^^h| Vw}, can be written as 

P„(t) = exp[-i(r!LBL + r!HBH)t] P^ exp[i(OLBL + rinBH)*]. (23) 

Because Pu> and commute, they can be simultaneously diagonalized by a unitary matrix X. Therefore, we have 
XHi^yC = diag[ft,(^,i, h^^2, ft-w.a], where h^^i < huj,2 < h^,z, and XP^^X''' = diag[pw, 1,^,^,2, — +^^,2)]- When /i varies 
slowly with t, the adiabatic precession solution can be obtained from the static precession solutions (with different /i) 
by using the adiabatic ansatz [82] 

In the two-flavor mixing scenario, the adiabatic ansatz corresponds to the assumption that P^^ remains parallel to 
H^. If neutrinos follows the adiabatic precession solution, then as ^ — > 0, P^j becomes a diagonal matrix in the mass 
basis (see Equations (22)). The diagonal elements of Pcj|^=o = PwU=o are p^,!, Puj.2 and —{puj.i +^^^,2), and these 
elements have the same order of appearance as do the diagonal elements of Ht^|^=o = ('^l — + (^h — f^H)BH- 

For example, in the mass basis we will have Pcj|^=o = diag[p^^2,Pcja, +^^^,2)] if H^^l^^o = diag[/i^_2, /icj,i, ^^^,3], 

where h^^i < h^^2 < hi^,3- The fact that Pi^|^=o is diagonal in the mass basis implies that there can be multiple 
spectral swaps/splits in the final neutrino energy spectra. Because a = <C 1, these spectral swaps/splits form 

hierarchically and appear at different neutrino densities (see Figure 7). 
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FIG. 7: The evolution of the diagonal elements of the neutrino density matrices in the mass basis during the hierarchical 
formation of spectral swaps/splits in a neutrino gas. Initially Pt^, and the density matrix p^i, is the same across the whole 
energy spectrum u: £ [0, 2a;o] (the thin lines in the panels on the right). When jj,^ luq, all neutrinos are in the heaviest eigenstate 
(i.e., the state with the largest eigenvalue) of ~ /xD. For the NH case two spectral swaps/splits form hierarchically at ^ ~ ojo 
and n ~ aujQ, respectively. Curves in the left-hand panels give pui with various values of cj. At ijq 2> A* ^ '^^o the heaviest 
eigenstate of H ~ (cj — nH,o)BH + pD is jz^s) for neutrinos with uj > iln.o- Subsequently these neutrinos no longer participate 
in collective oscillations. At p <^ auuo the remaining neutrinos are in 11^2) or {vi), depending on whether cijl ~ auj is larger 
or smaller than fiL.o- In either case the state corresponds to the heaviest eigenstate of ~ (cjl — ^^l,o)Bl (ignoring the 
decoupled state li^s)). This leads to two spectral swaps/splits in the final neutrino energy spectra (the thick lines in the NH 
column in the right-hand panels). In the IH case, initially the P33 are set to for all lj. As a result l^a) is completely decoupled 
from the EoM and no spectral swap/split forms at /i ~ ojo. The dotted lines in the right-hand panels are computed by using 
the constancy of Tr(D • A3) and Tr(D • Ag) and by assuming that the spectral swaps/splits are infinitely sharp. Reprinted 
figures with permission from B. Dasgupta et al, Phys. Rev. D, Vol. 77, 113007, 2008 (Reference [83]). Copyright 2008 by the 
American Physical Society. 

V. ANISOTROPIC AND/OR INHOMOGENEOUS NEUTRINO GASES 

A. Kinematic decoherence of collective neutrino oscillations 

The oscillations of neutrinos with different momenta can become out of phase (i.e., the breakdown of collective 
oscillations). This has sometimes been called "kinematic decoherence". Of course, this is not to be confused with 
quantum decoherence, which can be induced by any neutrino scattering process that changes neutrino momentum. In 
a homogeneous and isotropic neutrino gas the condition for the kinematic decoherence is Aw 3> /i. In this limit the 
coupling among "spins" is not strong enough to maintain a collective motion and will precess about B with vacuum 
oscillation frequency lj. In Section II C we have seen that, when Alj ^, synchronized neutrino oscillations do not 
decohere kinematically, and \D\ is approximately constant because of energy conservation. Numerical simulations 
suggest that bipolar neutrino oscillations are also stable [53, 55], although as yet we know of no conservation law 
which could explain this. 

If collective neutrino oscillations also exist in an anisotropic environment, then the wave fronts of the oscillation 
waves of the neutrinos must coincide with one another. In the neutrino bulb model these wavefronts are spheres that 
co-center with the PNS. However, in this model a neutrino propagating along a non-radial trajectory would travel a 
distance longer than that of a radially propagating neutrino [48]. In other words, a non-radially propagating neutrino 
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appears to have a larger oscillation frequency lo' = uj/ cos'd along r (see Equation (12)). This will enlarge Aw and 
would require stronger neutrino self-coupling to maintain collective oscillations among neutrinos propagating along 
different trajectories. Note that the "magnetic field" Hi/u,^(j') generated by other "spins" is trajectory dependent. 
This adds another potential source of kinematic dccoherence. 

Reference [70] pointed out that a neutrino gas with an initially symmetric bipolar configuration (i.e., P^, = ^Pv) 
can experience quick kinematic decoherence even in the presence of a small anisotropy, and in this case both \Pu\ and 
|Pp| evolve towards 0. Later it was found that an asymmetric bipolar neutrino system (i.e., \Py\ = (1 + £)|^p| with 
e 7^ 0) may nor may not experience kinematic decoherence depending on the value of e [73]. For typical choices of 
other parameters in a neutrino bulb model, kinematic decoherence is suppressed if e ^ 0.3. 

In the above discussion we have ignored ordinary matter. In the presence of a large matter density, P^ tends to 
precess around —L with frequency lo' ~ — A + iwcos20v In the isotropic environment lo' — > wcos20v in the corotating 
frame that rotates about L with frequency A. Therefore, in this case when neutrinos experience collective oscillations 
ordinary matter can be ignored. In an anisotropic environment, such as the neutrino bulb model, we will have 
frequency uj' = (—A + a;cos26'v)/ cosi? along r. Clearly, it is not possible to remove the matter effect for all neutrino 
trajectories. This implies that collective neutrino oscillations would not exist in a region where the net number density 
of electrons is much larger than that of neutrinos [89] . 



B. Precession mode in the anisotropic environment 

In Section III we have seen that with e ^ the flavor pendulum possesses an internal spin. The existence of this 
internal spin makes it possible for the flavor pendulum to experience simultaneously a precession motion and a nutation 
motion. These two kinds of motion of the flavor pendulum correspond to the precession mode and the bipolar mode of 
neutrino oscillations. The findings in References [70] and [73] suggest that bipolar neutrino oscillations, which arc the 
most prominent when e is small, becomes non-collective in the anisotropic environment. These findings also suggest 
that the precession mode, which becomes important for cases with sufficiently large e, can remain collective in the 
anisotropic environment. This is partially confirmed by the fact that neutrino oscillations calculated in the single- angle 
scheme and in the multi-angle scheme possess common features, such as the spectral swap/split phenomenon. This 
phenomenon, as explained with the single-angle supernova scheme or in homogeneous and isotropic environments, 
results from neutrino oscillations in the collective precession mode. 

The EoM of the polarization vectors in a stationary environment can be written as [111, 113, 122] 

p • V/i,p(x) = [loB + i/,,,p(x)] X /L,p(x), (25) 

where we have assumed A = for now, all the polarization vectors are normalized by n,y(xo), the number density of 
neutrino species v at location xq, and 

i?,,,p(x) = y2GFn.(xo) / Au' I d2p'(l - p . p')fl',p'(x). (26) 

Like the isotropic-neutrino-gas case. Equation (25) also exhibits rotational symmetry about B. If {Pi^^p(x)| Vw, p} 
solves the EoM, then \P'^ ^(x)| Vw, p} also solves the EoM, where P'^ ^(x) is obtained from Ptj.p(x) by rotation about 

B by an arbitrary angle </>, and where is independent of w, p, and x. This symmetry leads to the conservation law 

[911 



/ dw /"d2p(/t.p-S)p 

J — OO J 



= 0. (27) 



The rotational symmetry of the EoM about B can lead to a collective precession mode for neutrino oscillations, as 
shown in Figure 8, even in an anisotropic environment. 

Let K be the wave vector of the collective neutrino oscillation wave in a stationary environment. If -Plj,p experiences 
pure precession, then this vector must precess about B with frequency p • K (in flavor space) as the corresponding 
neutrino propagates along its word line (in coordinate space). Similar to the isotropic-neutrino-gas case, this means 

that Puj^f, is parallel to = (cj — p • K)_B -I- H^^ p. Therefore, the corresponding neutrino or antineutrino must be 
in a mass eigenstate if [91] 

[w-p-K] » li?,,,pl. (28) 
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FIG. 8; The collective precession of /L,^ in a neutrino bulb model with ^ 1 and the inverted mass hierarchy (Am^ < 0). 
The bottom three panels show that for most neutrinos the "spin components" si oc ■ e^^' and S2 oc P[^,,j ■ 62^' oscillate 
in phase right after collective neutrino oscillations begin (for this case, at r ~ 40 km). The in-phase oscillations of si and S2 
and the constant relative phase of these components imply the collective precession of the polarization vectors. The neutrinos 
or antineutrinos that drop out of the collective precession mode at smaller radii are those with vacuum oscillation frequencies 
farther from the collective oscillation frequency and those propagating along trajectories with larger cosi?_r. The top panel 
shows energy- averaged components, (si) and (S2). This panel shows that the collective precession mode stands out when 
non-collective neutrino oscillations average to zero. Figure adapted from Figure 1 in Reference [91]. Copyright 2009 by the 
Institute of Physics. 



Equation (28) gives the criterion for when neutrinos or antineutrinos drop out of the collective precession mode and 
begin to oscillate incoherently with respect to other neutrinos. The neutrinos or antineutrinos which drop out in this 
way have oscillation frequencies u so different from the collective oscillation frequency that neutrino self-interaction 
is not strong enough to maintain the corresponding polarization vectors in collective precession. 

With the replacements uj ((u;cos26'v — A) and B — >■ —L the above discussion also applies to the case with a large 
A. With these replacements. Equation (28) indicates that collective precession is indeed suppressed when the number 
density of electrons is much larger than that of neutrinos. 
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FIG. 9: Schematic plot of the neutrino oscillation regimes in the core collapse supernova environment in the two-flavor mixing 
scenario with |Am^| ~ Am^tm- Near the PNS, and in the "synchronized regime" {r < R~^n), collective neutrino oscillations are 
suppressed by either the large matter density or the large neutrino fluxes themselves. In this regime neutrinos can still experience 
i/-enhanced MSW flavor transformation. Far away from the PNS (r > -R^n), where the neutrino fluxes are negligible, neutrinos 
experience either vacuum oscillations or conventional MSW flavor transformation, depending on the matter density and the 
energy of the neutrino. If R^^n > R^oiu then there exists a window (-R^n ^ r < RtoW ^^e "bipolar regime") where neutrinos 
experience collective neutrino oscillations and where spectral swaps/splits develop. The curves show the energy-averaged 
neutrino survival probabilities for electron neutrinos and antineutrinos in single-angle calculations chosen to be representative 
of some late-time supernova conditions. 



VI. COLLECTIVE NEUTRINO OSCILLATIONS IN SUPERNOVAE 



A. Neutrino oscillation regimes 



With the picture for neutrino oscillations developed above, we can utilize the strength of neutrino self-interaction 
to sketch out neutrino oscillation regimes. These are shown in Figure 9. We will focus the following discussion on the 
two flavor mixing scenario with |Ato^| ~ Anil^^^ and 9^ ~ Oi^. We designate -R^n as the radius closest to the PNS 
where collective oscillations set in. Likewise, R^^n is the outer radius where collective oscillations cease. In both the 



NH and IH cases, can be estimated from the condition 



V2GFr...(i?+„)^^^, (29) 



where ni>_. (r) is the total number density of the antineutrinos at r that are initially Vg at the neutrino sphere, and 
is the average energy of these antineutrinos. Here we have chosen to use Vf, as the representative neutrino 

species, and we have used (^E^^f estimate of the frequency spread of the neutrino spectrum. 

In the NH case, assuming a fully-synchronized neutrino system, is approximately where [63] 

^/2GFr^e(i?7„ll) ~ ^ync (30) 

(In the large neutrino flux limit, the behavior of a neutrino system experiencing i^-enhanced MSW flavor transformation 
is the same as that of a synchronized neutrino system [74].) We note that for an iron core collapse supernova at early 
times and in the NH case, R^^n < -R^n, meaning that the collective neutrino oscillation regime does not exist. At later 
times in these models a collective neutrino oscillation regime can appear when the matter density becomes relatively 
small. In contrast, collective neutrino oscillations can occur at early epochs in an 0-Ne-Mg core collapse supernova. 
This is because the progenitors of these supernovae have relatively lower masses (8-I2M0), and therefore, after core 
bounce, the PNS has a dilute, lower density envelope [123, 124]. 

In the IH case the value of iJ^n is the larger of the values extracted from the following two conditions: 

1 Am^ 

"-(^-)=^(7TTi^71^ and n,.(i?,-,)...(i.,-,). (31) 
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The first condition in Equation (31) is based on the flavor pendulum model (See Equation (17)) and the second 
condition is in accord with the discussion in Section VA. In the very early epochs of an iron core collapse supernova, 
the collective neutrino oscillation regime may not exist in the IH case because of the presence of a large matter density. 

B. Effects of collective neutrino oscillations 

The supernova neutrino energy spectra can be dramatically modified by collective neutrino oscillations. Perhaps 
the most prominent feature to arise from collective neutrino oscillations is the spectral swap/split. For example, see 
Figure 1 and its caption. Obtaining this feature can depend on supernova conditions. For example, the calculations 
used to produce Figure 1 would give no spectral swap/split if the matter density were too large [76, 89]. 

When Vf. and are the most abundant neutrino species at the neutrino sphere, some qualitative aspects of collective 
neutrino oscillations in supernovae can be understood using the flavor pendulum model. In the IH case, the initial 
configuration of the flavor pendulum is near its highest position and is unstable when the neutrino flux is below 
some critical value. (See the discussion in Section IV B.) This implies that collective neutrino oscillations and their 
effects on neutrino energy spectra are relatively insensitive to the matter density or the exact value of ^13. In the 
NH case, however, the flavor pendulum would be near its stable configuration unless MSW flavor transformation 
displaces it from this position. As a result, collective neutrino oscillations and their effects depend on the efficiency 
of the {/-enhanced MSW flavor transformation. In turn, the efficiency of this flavor transformation is sensitive to 
both the matter profile and the value of ^13. This general picture is confirmed by multi-angle calculations [75]. 
The results shown in Reference [75] suggest that it may be possible to resolve the neutrino mass hierarchy with an 
observed supernova neutrino signal, even if the absolute neutrino masses and/or ^13 are too small to be measured in 
the laboratory. 

Neutrino signals detected at very early times after core bounce may be important probes of neutrino mixing. For 
iron core collapse supernovae, the neutrino spectra at the very early times, are modified by the conventional MSW 
mechanism and can be easily calculated [125]. For 0-Ne-Mg core collapse supernovae, however, the matter density is 
so low above the PNS that even at these very early times, where the v^, luminosity can be very large (~ 10^'^ erg/s), 
collective neutrino oscillations can create step-like features (swaps/splits) in the observed neutrino energy spectra 
[78]. If these features are observed in a Galactic supernova neutrino burst, then they could serve as diagnostics of the 
neutrino mass hierarchy and, in the NH case, even provide a measure of ^13. 

However, wc note that the spectral swap/split phenomenon is sensitive to the neutrino luminosities and energy 
spectra at the neutrino sphere. Presently there arc rather large uncertainties in these quantities, especially for the 
late-time supernova environment. We also note that very large neutrino fluxes do not necessarily imply a large 
neutrino oscillation effect. In fact, in the IH case a larger neutrino luminosity pushes -R^u to a larger radius. 

Of course, collective neutrino oscillations are not the only way that neutrinos may experience flavor transformation in 
supernovae. As in the Sun, in supernovae neutrinos also can experience the conventional MSW flavor transformation. 
However, the matter profile in supernovae may not be smooth (e.g., because of shocks or turbulence) near MSW 
resonance regions. This can produce some interesting phenomena [99, 126-132]. 

In principle, the alteration of supernova neutrino energy spectra by collective neutrino oscillations and/or other 
processes of neutrino flavor transformation could affect supernova dynamics, shock reheating, and nucleosynthesis in 
neutrino-heated ejecta [39, 46, 48, 106, 133]. However, modeling of this physics in realistic supernova conditions is at 
present primitive. 

VII. SUMMARY AND OPEN ISSUES 

In this review we have discussed collective neutrino oscillations in a simple bipolar neutrino system, in homogeneous 
and isotropic neutrino gases, and in anisotropic neutrino gases. The simple bipolar neutrino system (described by the 
flavor pendulum) is the simplest of these and is solvable analytically. It can be used to understand many qualitative 
features of collective neutrino oscillations in supernovae. The single-angle scheme essentially treats supernova neutrinos 
as a homogeneous and isotropic gas, and adiabatic neutrino flavor transformation in such a gas can be used to 
understand the spectral swap/split phenomenon in the supernova environment. An anisotropic neutrino gas can 
possess unique characteristics (e.g., suppression of collective flavor oscillations by large matter density). This makes it 
an important target for study because realistic physical environments such as core collapse supernovae can be highly 
anisotropic. 

We have covered some of the basic properties of collective neutrino oscillations, with emphasis on the two-flavor 
mixing case. There is much about the collective neutrino oscillation phenomenon which remains to be understood. 



18 



Antineutrinos Neutrinos 
— I 1 1 1 1 1 1 1 r 




10 20 30 40 10 20 30 40 50 
Energy [MeV] Energy [MeV] 



FIG. 10: Multiple spectral swaps/splits in a two-flavor single-angle calculation with [Am'^l ~ Am^tm and 6^, <^ 1. The solid and 
dashed curves are the initial and final energy spectra, respectively. The red and blue curves are for e and /i fiavors, respectively. 
The shaded region mark the energy ranges where spectral swaps occur. Unlike previous numerical calculations, this calculation 
assumes that i/^ and v^, rather than and u^, are the most abundant neutrino species at the neutrino sphere. Reprinted 
figure with permission from B. Dasgupta et al, Phys. Rev. Lett., Vol. 103, 051105, 2009 (Reference [101]). Copyright 2009 by 
the American Physical Society. 

For example, in Section IV A we skipped over the adiabaticity condition when we discussed adiabatic neutrino flavor 
transformation. Adiabaticity criteria developed so far [77] are difhcult to use in practice. 

There are other open issues in our current understanding of collective neutrino oscillations in supernovae. For 
example, it was recently shown that there can exist multiple spectral swaps/splits in the final neutrino energy spectra 
[76, 101, 103] (Figure 10). The appearance of these features depend on, however, the luminosities and energy spectra of 
the different neutrino species at the neutrino sphere. This finding cannot be explained by a grand collective precession 
mode in which all neutrinos and antineutrinos participate. In an even more recent paper [108] it was reported that 
qualitatively different results could appear in some collective neutrino oscillation scenarios depending on whether two- 
flavor mixing or full three-flavor mixing is used. All of these discoveries point up the need for a systematic study of 
supernova neutrino oscillations with neutrino self-interaction and full three-flavor mixing at various supernova epochs 
where neutrino luminosities and energy spectra can be different. In addition, numerical simulations of supernova 
explosions have shown clearly that realistic supernova environments may be highly anisotropic and inhomogeneous 
[28, 30-35, 134]. It would be interesting to see how density inhomogeneities arising from, e.g., shock waves and 
turbulence might affect collective neutrino oscillations. It remains a towering numerical challenge to integrate fully 
hydrodynamic, three-dimensional supernova models with calculations of neutrino oscillations that include neutrino 
self-interaction. However, such self-consistent integration may be necessary in order to study the interplay between 
supernova physics and neutrino oscillations. In particular, a full understanding of neutrino oscillations in supernovae 
may hold the key to deciphering the neutrino signal from a future Galactic supernova. The stakes are high. Deciphering 
a supernova neutrino signal could provide important insights into supernova astrophysics. It also could provide key 
insights into fundamental neutrino properties and these insights could be complementary to those sought by future 
neutrino experiments. 
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